Violation of finite-size scaling for the free energy near criticality 
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The singular part of the finite-size free energy density f s of the O(n) symmetric p 4 field theory is 
calculated for confined geometries of linear size L with periodic boundary conditions in the large-n 
limit and with Dirichlet boundary conditions in one-loop order. We find that both a sharp cutoff 
and a subleading long-range interaction cause a non-universal L dependence of f 3 near T c . For 
film geometry this implies a non-universal critical Casimir force with an algebraic L dependence 
that dominates the exponential finite-size scaling behavior above T c for both periodic and Dirichlet 
boundary conditions. 
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The concept of universal finite-size scaling has played 
an important role in the investigation of finite-size effects 
near critical points over the last decades . Consider 
the free-energy density f(t, L) of a finite system at the 
reduced temperature t = (T — T c ) /T c > and at van- 
ishing external field in a c?-dimensional cubic geometry 
of volume L d with periodic boundary conditions (pbc). 
It is well known that, for small t, the bulk free energy 
density ft = f(t, oo) can be decomposed as 



fb(t) = f bs (t) + f (t) 



(1) 



where fbs{t) denotes the singular part of fj, and where 
the regular part fo(t) can be identified unambiguously. 
According to Privman and Fisher ^,|| the singular part 
of the finite-size free-energy density may be defined by 



f.(t,L) = f(t,L)-f (t) 



(2) 



where fo is independent of L. The finite-size scaling 
hypothesis asserts that, below the upper critical dimen- 
sion d = 4 and in the absence of long-range interactions, 
f s (t,L) has the structure |J-||] 



f s (t,L) = L- d T{L/Z) 



(3) 



where J~(x) is a universal scaling function and £(t) is 
the bulk correlation length. Both £ and L are assumed 
to be sufficiently large compared to microscopic lengths 
(for example, the lattice spacing a of lattice models, the 
inverse cutoff A -1 of field theories, or the length scale 
of subleading long-range interactions). Eq. (||) includes 
the bulk limit f 3 (t, oo) — fp s (t) = Y£, d with a universal 
amplitude Y. Eqs. ([[]) - ([|) are expected to remain valid 
also for non-cubic geometries where the scaling function 
T(x) depends on the geometry and on the universality 
class of the bulk critical point but not on 5, or A and not 
on other interaction details 0]. In particular sublead- 
ing long-range interactions (such as van der Waals forces 
in fluids) that do not affect the universal bulk critical 



behavior are assumed to yield only negligible finite-size 
corrections to T(x). 

As a consequence, universal finite-size scaling proper- 
ties are generally believed to hold for observable quan- 
tities derived from f s {t,L), such as the critical Casimir 
force Fc in film geometry 

F c = -dr x (t,L)/dL (4) 

where the excess free energy per unit area is given by 

f ex (t,L)=Lf(t,L)-Lf b (t). (5) 

Near bulk criticality Eqs. (||) - @ yield 

F c ^,L) = L- d X c (L/d) (6) 

where 



X c {x) = (d- l)T(x) - x T'(x) + Yx d 



(7) 



with T'{x) = dJ-(x) / dx. The universal scaling 
structure of Eqs. @ and (0) has been confirmed 
by renormalization-group (RG) and model calculations 
Jl2|— and is predicted to be valid also for Dirichlet 
boundary conditions (Dbc) C^] . 

In this Letter we show that the scaling predictions of 
Eqs. (||) and (||) are significantly less generally valid than 
anticipated previously. On the basis of exact results for 
the O(n) symmetric ip 4 field theory in the large-n limit 
and on the basis of one-loop results we shall identify 
two sources in the ip 4 Hamiltonian that cause non-scaling 
finite-size effects: (i) a short-range interaction term ~ k 2 
with a sharp cutoff A in k space, (ii) an additional sub- 
leading long-range interaction term ~ 6|k| <7 ,2 < a < 4. 
In the latter case we consider both pbc and Dbc. Specif- 
ically, we find for 2 < d < 4 that Eqs. (||) and (||) must 
be complemented as 

f.(t,L,A) = L- 2 A d - 2 ^(r'A- 1 ) + L- d T(L/0 (8) 
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for the case (i), and 

F &L,b) = -bL- d+2 -° B(L/H) + L- d X c (L/0 (9) 

for the case (ii), respectively, where the function $ has 
a finite critical value $(0) > and where the function 
B(L/£) has a non- exponential decay ~ (L/£) -2 above T c 
for both pbc and Dbc. This implies (i) that the non- 
scaling L~ 2 term in Eq. (g) exhibits a dominant size de- 
pendence compared to the L~ d scaling term and (ii) that 
the non-universal term proportional to b in Eq. (M) im- 
plies an algebraic L dependence ~ bl; 2 L~ d ~ a that domi- 
nates the exponential finite-size scaling term of Xc above 
T c for both pbc and Dbc. By contrast, for the tp 4 lattice 
model with short-range interaction, we find that Eqs. (|^) 
and (||) are indeed valid except that for L 3> £ above T c 
the exponential scaling arguments of T and Xc must be 
formulated in terms of the lattice-dependent "exponen- 
tial" correlation length [^5| , ^6| . 

Non-negligible cutoff effects jl7],(l8| and non-universal 
finite-size effects due to subleading long-range interac- 
tions fl9| were already found previously for the finite- 
size susceptibility. These effects, however, were restricted 
to the regime L ^> £ close to the bulk limit above T c . 
The new non-scaling finite-size effect (i) exhibited in Eq. 
(g) is significantly more general in that it is pertinent to 
the entire — L^ 1 plane including the central finite- 
size regime £ 3> L. In particular, this effect exists at 
T c where / s (0, L,A) has a non- universal leading ampli- 
tude L~ 2 A d ~ 2 $(0) for d > 2. This implies the non- 
universality of the critical Casimir force L~ 2 A d ~ 2 $(0) 
in film geometry. 

Furthermore, the new non-scaling effect (ii) exhibited 
in Eq. (^|) has relevant physical consequences in systems 
with subleading long-range interactions. These conse- 
quences are significantly more important than those con- 
sidered previously for the finite-size susceptibility for pbc 
fl7| |l9| . The latter are of limited physical relevance since 
in real systems they are dominated by surface terms of 
0(L _1 ). In this Letter we predict a leading non-scaling 
effect on Fc not only for model systems with pbc but also 
for real systems with Dbc. 

We start from the standard tp 4 continuum Hamiltonian 



H 



d d x 



ro¥> 2 + 7; (V<p) 2 + u a (<p 2 f 



(10) 



with ro = ro c + a$t for the n-component field y(x) m 
a partially confined L d x oo d ~ d geometry with periodic 
boundary conditions in d' < d dimensions. This model 
requires a specification of the x dependence of </?(x) at 
short distances. We decompose the vector x as (y, z) 
where z denotes the coordinates in the d 1 confined direc- 
tions. We consider two cases (i) and (ii). 

Case (i) : We assume pbc and a sharp cutoff A, 
i.e., we assume that the Fourier amplitudes <^ p . q of 
(ys(y,z) = E p L £p,q e^ pz+qy ) are restricted to 

wave vectors p and q with components pj and qj in the 



range —A < pj < A and \qj\ < A. Here J q stands for 

(2ir)~ d+d f d d ~ d q , and ^ p runs over pj — 2nmj/L 
with m,j = 0, ±1,±2, .... The question can be raised 
whether there exists a non-negligible cutoff dependence 
of the finite-size free energy density per component (di- 
vided by fegT), 



fd,d'(t,L,A) = -n 



l L- d Jim L- d+d ' luZ, 

L — >oo 



d.d> 



(ii) 



where 



Z d , d >(t,L,L,A) = Y[ 

k,q 



d ^p, q 

A(d-2)/2 



cxp(-H) (12) 



is the dimensionless partition function of a L d x L d d 
geometry. For comparison we shall also consider the free 
energy density f(t, L, a) of the (p 4 lattice model 



H = a 



(13) 



with a nearest-neighbor coupling J on a simple-cubic lat- 
tice with a lattice spacing a. The factor (fc^T) -1 is ab- 
sorbed in H and H . 

We shall answer this question in the exactly solvable 
limit n — ► oo at fixed uqh where the free energy density 
is f§ 

fd,d> (t, L, A) = - h d \nn - {r ° ~ X ^ 
2 louon 

+ ^ d T/ ln [A- 2 (x- 1 + p 2 + q 2 )] ■ (14) 



Here x 1 is determined implicitly by 



X 1 = r Q + 4u Q n L d ' 



£ 



(15) 



The bulk free energy /& and bulk susceptibility \b above 
T c are obtained by the replacement L~ d J q ~* Jk' 
and the critical point is determined by ro = r$ c = 
—Auqu J k k~ 2 where k = (p,q) . The bulk correlation 

length above T c is £ = x^ 2 = £,ot~ v where v = (rf — 2) _1 . 
The regular part of fb reads /q = c\ A d — r 2 /(16uon) 
where c\ is a d dependent constant. The singular part of 
fb above T c is fb s = Y£, d with the universal amplitude 
Y = (d - 2) A d / [2d(4 - d)] where A d = 2 2 - d ^ d l 2 [d - 
2) _1 r(3 - d/2). For the singular part f s = f d ,d> - fo of 
the finite-size free energy above and at T c we find the 
form of Eq. @ with the leading non-scaling part 



6(2tt) 



d-2 



dy 



dq e 



d-i 



exp[-(l + r 2 A- 2 )2/] 



(16) 
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and the subleading universal scaling part 



2(4 - d) 



(L/H) d - 2 P 2 -~ d P d 



1 fdy f Fx 

2 J y \\y 

o 



d-d' 



W d ,{y)e- p y' A * (17) 
where P(L/£) is determined implicitly by 



P 



d-1 



(L/0 



d-2 



A-d 
An 2 A d 



dy 



x W d >(y)e 



o 

-P 2 y/4ir 2 



i-d' 



Wd(y)=Ln - E 



(18) 



(19) 



This result remains valid also for t < after replacing 
the terms (L/£) d - 2 in Eqs. @ and (Jl|) by t(L/£ ) d - 2 
and after dropping the term — ^ _2 A _2 ^in the exponent 
of Eq. ( Jl6[ ) . We have confirmed the structure of Eq. (|^) 
also for the <p 4 theory with finite n within a one-loop RG 
calculation at finite A which yields the same form of the 
function 3>d,d>(£,~ 1 A~ 1 ) as in Eq. (|l6|). Thus the univer- 
sal scaling form (||) is invalid for the tp 4 field theory with 
pbc and with a sharp cutoff, both for T > T c and for 
T<T C . 

These results have a significant consequence for the 
critical Casimir effect. Instead of Eq. (||) we obtain from 
Eqs. (|) and (jl|) - @ in film geometry (d' = 1) 

F c (t L, A) = L- 2 A- 2 $ d4 (r'A- 1 ) + L- d X c (L/£) . 

(20) 

Thus the critical Casimir force has a leading non- 
universal term ~ L~ 2 , in addition to the subleading uni- 
versal terms ~ L~ d of previous theories |9|-|l^], both for 
T > T c and for T < T c . 

We have also calculated fd.d' and Fc for the lattice 
Hamiltonian ( |l3| ) and for the continuum Hamiltonian 
( |i"o| ) with a smooth cutoff in the large-n limit. In both 
cases the scaling form (||) is found to be valid. For the 
lattice model, however, the second-moment bulk corre- 
lation length £ in the argument of T must be replaced 
by the lattice-dependent exponential correlation length 
pB|^6[ . Specifically we find, at fixed t > 0, the exponen- 
tial large-!/ behavior 

f s (t,L,a) - f bs = -d'(L/2^i) (rf - 1)/2 L- d exp(-L/6) 

(21) 

where £i = (a/2) [arcsinh(a/2^)] _1 is the exponential 
correlation length in the direction of one of the cubic 



axes. Note that the non-universal dependence of £i on a 
is non- negligible in the exponent of ( pTj ) [ [IS) . 

The sensitivity of f s (t,L,A) and Fc with respect to 
the cutoff procedure can be explained in terms of a cor- 
responding sensitivity of the bulk correlation function 
G(x) = < cp(x.)cp(0) > in the range |x| » £ @. For 
example, for the tp 4 continuum Hamiltonian fllpp with 
an isotropic sharp cutoff |k| < Awe find, in the large-n 
limit, the oscillatory power-law decay above T c 

d-2m_„^-(d+i)/2 sin [Ax-?r(d- l)/4] 



G(x) = 2 A d -\2nxA) 
+ o( e~ x ^ 



l + £- 2 A- 



(22) 



for large x = |x| 3> £ corresponding to the existence of 
long-range spatial correlations which dominate the expo- 
nential scaling dependence ~ e~ x ^. By contrast, G(x) 
has an exponential decay for the lattice model ([13]) with 
purely short-range interaction [fl5[ . An exponential de- 
cay of G(x) is also valid for the continuum model (JTc" 
with a smooth cutoff jl5| . 

The non-universal cutoff effects on f s , Fq and on G(x) 
described above are a consequence of the long-range cor- 
relations induced by the sharp-cutoff procedure in the 
presence of pbc. We consider these consequences not 
only as a mathematical artifact but rather as an impor- 
tant signal for a serious lack of universality in physical 
systems. We substantiate this interpretation by demon- 
strating that corresponding violations of finite-size scal- 
ing should indeed exist in physical systems with more 
realistic interactions and boundary conditions. 

Case (ii) : We assume the existence of a sublead- 
ing long-range interaction in the continuum tp Hamil- 
tonian H which in the Fourier representation has the 
form 6|k| a with 2 < a < 4, in addition to the short-range 
term k 2 . It is well known that the subleading interac- 
tion ~ |k| CT corresponds to a spatial interaction potential 



V(x) 



-d—a 



that does not change the universal bulk 



critical behavior |2lj]. Interactions of this type exist in 
real fluids. As pointed out by Dantchev and Rudnick 
M|, the presence of this interaction yields leading non- 
scaling finite-size effects on the susceptibility x f° r the 
case of pbc in the regime L £ above T c , similar to 
those found for a sharp cutoff |L7]jll|. 

In real systems with non-periodic boundary conditions, 
however, these non-scaling finite-size effects become only 
subleading corrections that are dominated by the surface 
terms of x of 0(L^ 4 ). In the following we show that the 
situation is fundamentally different for Fq which, by def- 
inition, does not contain contributions of 0{L~ l ) arising 
from the 0(L _1 ) part / of the free energy density. 

We consider film geometry and first assume Dbc in 
the z direction corresponding to <p(y,0) = ip(y,L) = 
0, i.e., we assume that E p in the Fourier representa- 
tion of tp(y,z) = L^ 1 J2 P Jq 0p,q e jqy sin(pz) runs over 
p = mn/L,m = 1,2,.... Such boundary conditions are 
relevant for the superfluid transition of 4 He |23]. The 
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presence of the subleading interaction £>|k|' T implies, for 
L > ^ above T c , a non- universal term ~ b in 

f,(t,L,b)-f = -bL- d+2 -^(L/0 + L- d G(L/0 

(23) 

where Q(L/^) is the known p2[ universal scaling function 
for purely short-range interaction. The scaling function 
Q{L/^) has an exponential large-L behavior ]l2f . By con- 
trast we find that \&(L/£) has an algebraic L dependence. 
Performing a one-loop RG calculation we obtain 

oc 

*(L/0 = \ (27r) CT - 4 J dx(l + x-^)*(x), (24) 

{L/ay- 

oo d—1 

*M = | * » <2 -'" 2 e-»' 4 "* (J?) x 

where 7*(z, x) = £~ z J X die - * t 2 ' 1 / J °° dte~* is 
the incomplete Gamma function and 

— -.00 

I - 2 £ cx K-H- (26) 

n— — 00 

We have found that cutoff effects are negligible for the 
function i&(L/£). At fixed £, the large-L behavior is 
*(L/£) ~ ( L /0~ 2 - Eq. © yields the following form 

S(L/0 = (d-3 + ff )f(L/f)-(i/p'(I/0 (27) 

for the non-universal contribution to H7 in Eq. @ . The 
crucial consequence is that the leading critical tempera- 
ture dependence ~ 6£ 2 L~ d ~ CT of Fc for L > £ above T c 
is algebraic and non-universal whereas the critical tem- 
perature dependence of the scaling part Xc{L/£) derived 
from C?(L/£) is exponential and universal Jl^ ]. This may 
have significant consequences for the interpretation of ex- 
isting |23| and future experimental data in fluids. 

For comparison we finally present our result for 
for film geometry in the presence of pbc. In one-loop 
order we obtain for ^ p i, c (L/^) the same form as given 

for in E q s - @ and © but with Wi(y) re- 

placed by W\(y), Eq. (JTg|) . For the large-L behavior 
we find if?pb c (L/£) ~ (L/£;)~ 2 which dominates the ex- 
ponential scaling dependence of Ap. Our prediction of 
a non-exponential non-scaling effect on Fc above T c for 
pbc can be tested by Monte Carlo simulations (lO) . 

The next stage of the theory would entail a quanti- 
tative determination of the non-universal interaction pa- 
rameter b for a specific system. This is, of course, beyond 
the scope of the present paper. 
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